T 


AD-A132  017 


UNCLASSIFIED 


PRELIMINARV  DATA  PROCESSING  PLAN  FOR  THE  THERMAL  PLASMA 
EXPERIMENT  ON  THE  HILAT  SATELLITE(U)  AIR  FORCE 
GEOPHVSICS  LAB  HANSCOM  AFB  Mf)  FJ  RICH  ET  AL. 

21  MAR  82  AFGL-TR-83-0091  F/G  5/2 


1/i 


NL 


OTIC  FILE  COPY 


ff 


AFGL-TR-83-0091 

INSTRUMENTATION  PAPERS,  NO.  317 


[S 

A 


V*  1 


Preliminary  Data  Processing  Plan  for  the  Thermal 
Plasma  Experiment  on  the  HILAT  Satellite 

F.  J.  RICH 
R.  A.  HEELIS 


31  March  1983 


Approved  (or  public  release;  distribution  unlimited. 


SPACE  PHYSICS  DIVISION  PROJECT  7601 

AIR  FORCE  GEOPHYSICS  LABORATORY 


HANSCOM  AFB,  MASSACHUSETTS  01731 

AIR  FORCE  SYSTEMS  COMMAND,  USAF 


83 


i 


DR.  ALVA  T.  STAIR,  Jr^ 

Chief  Scientist 

Qualified  requestors  may  obtain  additional  copies  from  the 
Defense  Technical  Information  Center.  All  others  should  apply 
to  the  National  Technical  Information  Service. 


EDITION  OP  '  NOV  65  IS  OBSOLETE 


_  Unclassified 

SECURITY  CLASSiFi  at, on  OP  This  PAGE  »>irn  D»tm  Entered) 


_ Unclassified _ 

SECURITY  CLASSIFICATION  QP  THIS  PAGECW?>fi  Data  Ent»r»d) 

20.  Abstract  (Contd) 

information.  Refinements  to  the  algorithm  will  be  required  before  the  final 
results  are  released. 


Contents 


1.  GENERAL  INFORMATION 


2.  DATA  FORMAT  8 

2.  1  Telemetered  (TM)  Data  C 

2.2  DataStorage/Readout  t> 

2.3  SUBCOM  Word  10 

2.4  Timing  10 

3.  DRIFTMETER  DATA  PROCESSING  11 

3.  1  Theory  of  Operation  n 

3.2  Operation  of  Driftmeter  14 

3.2.1  Re-Zero  Measurement  15 

3.2.2  Offset  Measurement  111 

3.2.3  Difference  Measurement  17 

3.2.4  TM  Outputs  17 

3.3  Data  Analysis  13 

4.  RETARDING  POTENTIAL  ANALYZER  (R PA)  DATA 

PROCESSING  21 

4.  1  Theory  of  Operation  21 

4.2  Operation  of  RPA  2  5 

4.3  Data  Analysis  30 

5.  ELECTRON  SENSOR  DATA  PROCESSING  3! 

5.  1  Theory  of  Operation  31 

5.2  Operation  of  Electron  Sensor  5  5 

5.3  Data  Analysis  3  5 

5.3.1  Density  Mode  Analysis  3.j 

5.3.2  Swept  Applied  Potential  Mode  Analysis  39 

II.  CONCLUSIONS  40 


LIST  OF  SYMBOLS 


43 


Illustrations 


1.  DNA/HILAT  Configurations  0 

2.  Experiment  Block  Diagram  7 

3.  HILAT  Driftmeter  Collector  and  Aperture  Configuration  12 

4.  Cross-section  of  Retarding  Potential  Analyzer  for 

DNA/HILAT  Satellite  22 

5.  An  Example  of  Ion  RPA  and  Spherical  Electron  Sensor 

Data  as  a  Function  of  Applied  Voltage  at  840-km 

Altitude  24 

6.  Schematic  of  Applied  Voltage  Sequence  for  RPA  and 

Electron  Sensor  2(i 

7.  Example  of  log^Q<-Ie)  vs  <£p  for  HILAT  Electron  Sensor  34 

8.  DNA/HILAT  Langmuir  Probe  Filter  Bank  Response 

Functions  3  7 


Tables 


1.  Data  Sampling  Format  for  Frames  0  to  123  Thermal 

Plasma  Experiment  —  DNA/HILAT  Satellite  9 

2.  Sequence  of  Values  in  Driftmeter  Measurement 

Word  (I DM)  19 

3.  Applied  Voltages  for  DNA/HILAT  27 

4.  Calibration  of  RPA  Log  Electrometer  29 

5.  Calibration  of  Electron  Sensor  Log  Electrometer  3(i 

6.  Filtered  Output  vs  M- Ie)/(-Ip)  38 

7.  Summary  of  Conversions  of  Raw  Data  From  TM  Units  41 

8.  Summary  of  Preliminary  Calculated  Results  41 


4 


Preliminary  Data  Processing  Plan  for  the 
Thermal  Plasma  Experiment  on  the  HILAT  Satellite 


1.  GENERAL  INFORMATION 

The  Thermal  Plasma  Experiment  is  one  of  five  experiments  to  be  flown  on  the 
HILAT  Satellite  sponsored  by  the  Defense  Nuclear  Agency.  The  HILAT  spacecraft 
is  scheduled  to  be  launched  into  a  830-km  circular  orbit  with  an  inclination  if 
82.2°.  The  +x  axis  will  be  along  the  velocity  vector  (±10°).  The  -  "  ‘v 
along  the  radial  vector  (±10°).  The  prime  experiment  of  the  mission  is  the  r  ho 
beacon,  or  Wide  Band,  which  transmits  a  phase  coherent  signal  at  YHF,  lTHF, 
and  L-band  frequencies.  The  signal  is  used  to  probe  the  ionosphere  between  the 
satellite  and  ground  receivers  for  regions  that  cause  scintillation. 

The  other  experiments  onboard  are:  the  Thermal  Plasma  Experiment,  which 
is  described  in  detail  in  this  report;  the  Particle  Spectrometer,  which  measures 
electron  fluxes  from  20  eV  to  20  keV;  the  Fluxgate  Magnetometer,  which  is  used 
both  for  attitude  determination  and  for  detecting  field -aligned  currents;  and  the 
Auroral  Imager/Mapper,  which  images  the  ionosphere  below  the  spacecraft  in  the 
vacuum  ultraviolet  frequencies. 

The  Thermal  Plasma  Experiment  consists  of  three  sensors,  a  ground  plane, 
a  mounting  bracket  (to  hold  the  sensors  and  ground  plane),  and  an  electronics  box. 
The  configuration  of  the  experiment  on  the  spacecraft  is  shown  in  Figure  1.  The 
three  sensors  of  the  experiment  are  the  Retarding  Potential  Analyzer  (HPA),  the 

(Received  for  publication  25  March  1083) 


5 


Driftmeter  (DM),  and  the  electron  sensor.  The  normal  to  the  HPA  and  DM  aper¬ 
tures  are  parallel  to  the  spacecraft  +X  axis.  The  electron  sensor  boom  is  parallel 
to  the  spacecraft  -Y  axis.  The  ground  plane  is  parallel  to  the  spacecraft  Y-Z 
plane.  The  functional  block  diagram  of  the  experiment  is  shown  in  Figure  2. 


[  “  -  ~  . .  I 

Cii  I  -  WIDE'  bAND  AMTENNA 


Figure  1.  D Y.A /III FAT  Configurations.  (-experiment  deck  and  penthouse  sections 
Seen  from  the  forward  (+x)  side 


2.  DATA  FORMAT 


2.1  Telemetered  (TM)  Data 

.  RPA  electrometer  (RPA) 

.  RPA  sweep  monitor  (SWP) 

.  Ion  Driftmeter  measurements  (IDM) 

.  Frame  Counter  (COUNT)  -  counts  TM  frames  in  modulus  128  (0  to  127) 

.  IDM  log  electrometer  outputs  (LLA  and  LLB) 

.  Electron  probe  electrometer  (ELEC) 

.  Temperature  monitor  (TEMP) 

.  Sensor  potential  monitor  (SENPOT) 

.  Four  bandpass  filter  outputs  (FI,  F2,  F3,  F4) 

2.2  Data  Storage/Readout 

Satellite  operations  are  divided  into  TM  frames  of  1/2  sec  of  spacecraft  time. 
(1.  000  sec  in  spacecraft  time  =  0.  99950  sec  real  time.  All  time  measurements 
given  herein  are  in  spacecraft  time.)  The  electronics  for  the  Thermal  Plasma 
Experiment  starts  its  frame  with  the  start  of  the  read-out  gate  (ROG),  which  lags 
the  start -of-frame  (SOF)  pulse  by  375  msec.  The  experiment  frame  is  divided 
into  04  equal  timeslots  by  a  128-Hz  counter.  The  TM  words  are  sampled  accord¬ 
ing  to  the  timeslot  in  which  they  occur.  The  data  is  sampled  and  stored  5.  13  msec 
after  the  beginning  of  the  timeslot.  The  512  bits/frame  are  allocated  for  seconds 
0  to  01  (frames  0  to  123)  as  follows: 

Word  Content.  *  Words/Frame  H  Bits/Word  Total  #  Bits 


Table  1.  Data  Sampling  Format  for  Framed  0  to  123  Thermal  Plasma 
Experiment  —  DNA/HILAT  Satellite 


TM 

Word 

Count 

Time- 

slot 

Bits 

Contents 

TM 

Word 

Count 

Time¬ 

slot 

Bits 

Contents 

0 

0-7 

SUBCOM 

28 

32 

256-263 

SUBCOM 

1 

8-16 

ELEC 

29 

33 

264-272 

ELEC 

2 

17-26 

IDM 

30 

34 

273-282 

IDM 

3 

27-35 

RPA 

31 

35 

283-291 

RPA 

4 

3  6-44 

ELEC 

32 

37 

292-300 

ELEC 

5 

6 

45-54 

IDM 

33 

38 

301-310 

IDM 

6 

7 

55-63 

RPA 

34 

39 

311-319 

RPA 

7 

8 

64-71 

SWP 

35 

40 

320-327 

SWP 

8 

9 

72-80 

ELEC 

3  6 

41 

238-336 

EL  K( ' 

9 

10 

81-90 

IDM 

37 

42 

337-346 

IDM 

10 

11 

91-99 

RPA 

38 

43 

347-355 

RPA 

11 

13 

100-108 

ELEC 

39 

45 

356-364 

ELEC 

12 

14 

109-118 

IDM 

40 

46 

365-374 

I  DM 

13 

15 

119-127 

RPA 

41 

47 

375-383 

RPA 

14 

16 

128-135 

SUBCOM 

42 

48 

384-391 

SUBCOM 

15 

17 

136-144 

ELEC 

43 

49 

392-400 

ELEC 

16 

18 

145-154 

IDM 

44 

50 

401-410 

IDM 

17 

19 

155-163 

RPA 

45 

51 

411-419 

U  PA 

18 

164-172 

ELEC 

46 

53 

420-428 

ELEC 

19 

173-182 

IDM 

47 

54 

429-438 

IDM 

20 

183-191 

RPA 

48 

55 

439-447 

RPA 

21 

192-199 

SWP 

49 

56 

448-455 

SW  P 

22 

200-208 

ELEC 

50 

57 

456-464 

EL  EC 

23 

209-218 

IDM 

51 

58 

465-474 

IDM 

24 

27 

219-227 

RPA 

52 

59 

475-483 

RPA 

2  5 

MB 

228-236 

ELEC 

53 

61 

484-492 

ELEC 

20 

iKi;  -M 

237-246 

IDM 

54 

62 

493-502 

IDM 

27 

247-255 

RPA 

55 

63 

503 -5  i  i 

RPA 

The  512  bits/ frame*  are  allocated  for  seconds  62  and  63  (frames  124  to  127)  as 
follows: 

Word  Content  #  Words/Frame  #  Bits /Word  Total  »  Bits 

144 


I.  LB /L  LA 
ELEC 
SWP 
IDM 

SUBCOM 


10 

it; 

4 

p; 

_4 

56 


9 

9 

8 

10 

8 


144 
3  2 
1  40 
32 
512 


Total 


The  format  for  framed  124  to  127  are  the  same  as  frames  0  to  123,  except  that 
LLB  replaces  RPA  for  Til  words  3,  10,  17  etc.,  and  1.1. .4  replaces  I \ P A  for  Til 
words  6,  13,  20  etc. 

2.3  SUBCOM  Word 

The  data  in  the  SUBCOM  word  is  multiplexed  .as  follows: 


Word  Count 

Frame  « 

0,  4,  8, _ 124 

Frame  * 

1,  3,  5,  .  .  . 

,  127 

Frame  « 

2.  n,  10 . 120 

0 

COUNT 

11 

COUNT 

14 

I.  LA 

12 

LLB 

28 

ELEC 

13 

ELEC 

42 

TEMP 

F4 

SENPOT 

TEMP  (0  to  255  =  0.00  to  5.  10V)  is  the  output  of  the  temperature  monitor  inside 
the  electronics  package  for  the  experiment.  The  conversion  from  telemetry  to 
temperature  is: 

T  =  -40°C  +  TEMP  *  (0.020V/bit)*(25°C  V)  .  (1) 

SENPOT  (0  to  255  =  0.  00  to  5.  10V)  is  the  output  of  the  floating  bias  volt;  je  applied 
between  the  spacecraft  ground  and  the  experiment  ground.  The  conversion  from 
telemetry  to  potential  is: 

<iSFN  -  IB.  0V  -  SENPOT  «  4.0  V/V  .  (2) 

The  nominal  limits  are  "  0.  00V  (SENPOT  =  4.48V)  and  bs  =  14.0V 

(SENPOT  =  1.00V). 

2.4  Timing 

Two  data  storag.  registers  are  used  so  ll'.at  one  register  may  be  loaded  while 
the  other  is  read  out  by  the  TM.  Data  readout  will  lag  data  store.,  a  by  1  2  of  a 
spacecraft  second.  Each  word  is  serially  stored  as  it  occurs  in  the  data  format. 
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DRIFT  METER 


radially 

down 


Figure  3.  IIILAT  Driftmeter  Collector  and  Aperture  Configuration 


Likewise  the  current  to  plates  B  and  C  is 

IRC  =  K  *  (1/2  W  -  D  *  tan  (a))  .  (4) 

The  value  that  is  sent  to  the  Til  is  log  I^R  -  log  1^,^.  It  van  be  shown  that  for 
|  a  |  <:  2  5° 

1°glO(IAB/ICD)  *  k  *  tan  a  '  (f,) 

where 

k  *  -2.3  *  tan  a  *  W/(4  *  D)  . 

Then,  the  horizontal  component  of  the  plasma  drift  is 

V  --  (V1  +  V  ,)  *  tan  (a)  -v  •  V  ,  .  CO 

y  s  u  ~s 

Ideally,  will  be  determined  from  the  RPA  data.  However,  for  preliminary 
processing,  we  will  set  this  value  to  zero. 

The  procedure  for  obtaining  the  vertical  component  of  the  plasma  drift  is 
identical  except  the  TM  value  is  log  -  log  1CR.  Then,  using  Lq.  (5),  the 
vertical  component  of  the  plasma  drift  is 

V  =  (V1  +  V  ,)  *  tan  a  +  z  .  V  .  '7) 

z  s  d  ~  • — s 

We  can  estimate  the  component  of  the  drift  that  is  normal  to  the  sensor  aperture 
by  assuming  that  the  drift  parallel  to  the  magnetic  field  is  zero: 

V  .  -  -  (V  •  B  +  V  •  B  )/B  ,  (8) 

d  ~-7.  z  ~y  y  '  x  ’ 

where 

B  ,  B  ,  B  =  component  of  a  unit  vector  along  the  magnetic  field  vector, 
x  y  z 

A  series  of  screens  or  wire  grids  are  located  between  the  aperture  and  the 
collector  plates,  but  are  not  shown  in  Figure  3.  The  grids  have  electrostatic 
potentials  applied  to  them  to  perform  a  variety  of  functions.  These  include  the 
surpression  of  electron  currents  to  and  from  the  collector  plates,  exclusion  of 
light  ions  upon  occasion  as  explained  in  Section  3.2,  and  compension  of  the  ion 
trajectories  due  to  the  other  potentials.  As  a  result  of  the  grids,  the  effective 
value  of  D  and  the  actual  value  are  slightly  different. 
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3.2  Operation  of  Driftmeter 


There  are  two  log  electrometers  used  by  the  driftmeter.  Their  outputs  are 
referred  to  as  and  Vg,  which  are  placed  directly  into  the  TM  as  LLA  and  LLB. 
If  the  log  electrometers  ere  connected  to  plates  A  +  B  and  C  +  D,  the  outputs  are 

LLA  -  V1  -  Gl  *  log  IAfi  +  V1Q 

LLB  ='  V2  =  g2  *  !og  !cd  +  V20  -  <9> 

where 

Gj,  Gg  =  amplifier  gains  =  1.01 

Viq.  =  amplifier  zero  levels  s  10.62. 

We  assume  that  the  log  amplifiers  are  matched  so  that  G^  =  Gg  =  Gj  .  The  zero 
levels  are  set  so  that  V.  =  Vg  =  «  4.  56  V  when  I  =  10~9  A. 

Instead  of  sending  and  Vg  to  the  TM  for  use  in  determining  the  ion  drift 

velocities,  they  are  input  to  a  difference  electrometer.  This  output  is  amplified 
via  a  capacitor  with  a  voltage  Vp  across  it.  The  final  voltage  presented  to  the 
telemeter  is 

VD  =  G3  *  «V2  -  Vl>  +  VcJ  •  (10) 

The  value  of  is  set  once  each  4  sec  to  the  instantaneous  value  of  (V^  -  Vg),  so 
that  Vp  =  0  if  the  environment  is  changing  slowly.  Since  there  are  two  modes  of 
connecting  the  log  electrometers  when  measuring  the  ion  drifts,  there  are  two 
capacitor  voltages  stored;  one  for  the  horizontal  drifts  and  one  for  the  vertical 
drifts.  The  gain  of  the  difference  amplifier  has  four  levels  of  sensitivity  such 
that 


Gg  -  J  *  4(3"H)  ,  (1] 

where 

11  range  =  0,  1,  2,  3. 

The  difference  amplifier  output  is  translated  by  2.52  V  before  being  presented  to 
the  TM.  Thus,  the  difference  value,  D,  in  the  TM  equals  Vp.  +  2.52  V. 


The  values  of  R  and  D  are  combined  into  a  single  TM  word  so  ’hat 


D. 

i 

R . 

l 

value 

0 

TTr  M 

1  2  3  4  5  5  7 

l 

8  9 

bits 

M 

L 

C 

s 

B 

O 

B 

w  here 

i  =  index  within  a  frame  1,  2,  3,  ,  lf> 

Dj  -  difference  output  0  V  to  5.  10  V 

R^  =  I'ange  -  0,  1,  2,  3. 

The  driftmeter  operates  on  a  4-Sec  (eight-frame)  cycle.  In  addition  to  n'pi1.'- 
senting  the  ion  drift  velocity  measurements,  th--  first  eight  values  of  the  firs’  of 
;hese  eight  frames  represent  "re-zero"  and  "offset"  measurement;  which  ai  • 
described  in  the  next  section.  The  proper  frame  in  which  the  re-zero  and  oiV»<  t 
measurements  are  located  is  found  with  the  use  nr  the  f:  ante  counter,  COUNT, 

The  superscripts  R,  O.  and  D  used  in  Section  3.2.  1  denote  values  related  to 
re-zero,  offset,  and  drift  measurements. 

3.2.1  RE-ZERO  MEASUREMENTS 

The  re-zero  measurement  is  made  by  disconnecting  the  output  th,.  difference 
amplifier  from  the  T>.1  line  and  grounding  it  through  tin  capacitor  u  .  iui  m< 
horizontal  or  vertical  drift.  The  output  to  the  TM  during  a  re-zero  measurement 
is  the  translation  voltage  that  should  be  2.52  V.  Small  variations  from  the  nominal 
translation  is  anticipated  due  to  system  temperature  changes.  The  true  transla- 
tion  voltage,  D  ,  should  be  used  in  the  data  analysis.  The  main  purpose  of  th. 
re-zero  measurement  is  to  re-sot  the  level  of  V  used  in  Kq.  (10).  Using  tin 
example  of  Eq.  (9),  the  horizontal  capacitor  would  be  se:  io 

Vc  “  V1  '  V2  -  °L  *  +  (V10  -  V20>  •  (12) 

If  the  environment  is  unchanging,  the  following  drift  measurements  will  be  nearly 
equal  to  the  null  value  output  to  the  TM  during  the  re-zero  measurement.  Any 
changes  in  the  environment  will  be  measured  relative  to  the  value  of  V  set  during 
the  preceding  re-zero  measurement. 


15 


3.2.2  OFFSET  MEASUREMENT 


The  purpose  of  the  offset  measurement  is  to  find  the  voltage  placed  on  the 
capacitor  during  the  re-zero  measurement.  The  offset  measurement  is  made  by 
reversing  the  log  amplifier  input  to  the  difference  amplifier.  The  offset  measure¬ 
ment  is  made  immediately  following  a  re-zero  measurement  and  it  depends  upon 
the  assumption  that  the  ion  arrival  angle  and  the  outputs  of  the  log  amplifiers  do 
not  change  between  the  two  measurements.  The  result  is  that  the  output  of  the 
difference  amplifier,  using  the  example  of  Eq.  (9),  is 

VD°  =  G3  *  «Vl  -  v2)  +  vc] 

=  G3  *  [2  *  Gl  *  log(IAB/ICD>]  •  (13) 


By  substituting  Eqs.  (5)  and  (11)  into  Eq.  (13),  we  find  the  difference  output  value 
in  the  TM  during  a  re-zero  measurement  is 

D°  =  2  *  Gl  *  k  *  J  *  4<3'R  )  *  tan  a0  +  2.  52  V  .  (14) 

The  ion  drift  angle  during  the  offset  measurement  is  determined  from 


0  _  D°  -  2.  52  V _ 

tsn  3  -  T- 

2*GL*k*J*4  **(3  -  RU) 


(15) 


By  substituting  Eq.  (15)  into  Eqs.  (6)  or  (7),  the  ion  drift  velocity  component  is 
obtained. 

From  t lie  offset  measurement,  we  can  find  a  voltage  equal  to  the  voltagf  placed 
on  the  capacitor  during  the  re-zero  measurement  plus  any  mismatch  between  the 
log  amplifiers.  By  substituting  Eqs.  (13)  and  (14)  into  Eq.  (12),  we  obtain  the 
voltage  on  the  capacitor  as 


V 


D  -  2.  52  V 
2  *  G„° 


(V 


10 


'  V20' 


(16) 
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I  he  offset  value  to  be  used  in  the  following  difference  measurements  is 


0  =  Vc  +  (V20 


V  )  s  ■  P  -2-  52  V _ 

v  1  (V  0 

2  *  J  *  4**(3  -  R  ) 


3.2.3  DIFFERENCE  MEASUREMENT 

The  majority  of  values  entered  into  the  TM  for  I  DM  represent  tJ.- ■  ltff.-r. 
measurements.  For  the  differences  due  to  horizontal  ion  drifts,  plates  A  ;.tt  i  1 
are  connected  to  the  first  log  amplifier  and  C  and  D  are  connected  to  the  s-  ■ 

For  vertical  drifts,  plates  A  and  D  are  connected  to  the  first  log  amplifier  arid  B 
and  D  are  connected  to  the  second.  In  addition,  the  voltages  stored  on  the  cair  - 
itors  during  the  horizontal  and  vertical  re-zero  are  used  in  combination  with  tie- 
output  of  the  log  amplifiers  as  input  to  the  difference  amplifier.  From  Eqs.  (9), 
<101,  and  (11),  the  voltages  from  the  difference  amplifier  is 

V0°  =  (J  *  4**<3  -  R°))  *  (Vc  -  *  log(IAB/ICD)  +  (V20  -  V10»  .  (13) 

Using  the  previously -obtained  value  for  V  ,  Eq.  (IT),  the  difference  voltage  is 

VD°  =  (J  *  4**(3  -  R°))  *  (0  -  GL  *  log(IAB/ICD))  .  (19) 

Using  Eq.  (5),  the  difference  measurement  put  into  TM  is 

D.°  =  (J  *  4**(3  -  RtD))  *  (0  -  Gl  *  k  *  tan  a.D)  +  2.  52  ,  (20) 

where  the  subscript  i  is  used  to  denote  the  sample  within  the  frame.  The  drift 
angle  can  be  found  as  a  function  of  values  found  in  the  TM  by  rearranging  term:-'  in 
Eq.  (20) 

tan  a.D  =  (0  -  (DjD  -  2.  52)/(J  *  4**<3  -  R.D)))/(GL  *  k)  .  (21) 

3.2.4  TM  OUTPUTS 

As  previously  explained,  the  outputs  of  the  driftmeter  are  1DM  that  oceur 
10  times  per  frame  and  LEA  and  I, LB  that  are  in  the  SUBCOM  word.  The  operation 
of  the  driftmeter  represents  an  8-sec  cycle  composed  of  two  sub-cycles.  During 
the  first  4  sec  of  the  8-sec  cycle,  the  driftmeter  makes  32  horizontal  drift  meas¬ 
urements  per  second.  During  the  second  4  sec,  the  driftmeter  alternates  between 
horizontal  and  vertical  measurements,  making  16  measurements  per  second  of 
each. 
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The  first  eight  values  of  IDM  in  each  4-sec  period  represent  the  re-zero  and 
offset  measurements.  The  first  four  values  represent  re-zero  and  offset  meas¬ 
urements  made  with  a  series  of  electrostatic  potentials  on  the  screens  between 
the  aperture  and  the  collecting  plates.  These  potentials  are  designed  to  prevent 
the  H+  ions  from  reaching  the  collector  plates  without  affecting  the  trajectory  of 
the  0+  ions.  The  second  four  values  of  the  first  frame  of  data  are  re-zero  and 

-f*  "f* 

offset  measurements  made  with  a  potential  on  the  screen  set  to  pass  H  and  O 
ions  equally  without  affecting  their  trajectories  from  the  aperture  to  the  collecting 
plates.  The  screens  are  used  to  electrostatically  prevent  secondary  electrons 
from  leaving  the  collector  plate.  .All  drift  measurements  are  made  with  the  later 
set  of  potentials  on  the  screens. 

The  arrangement  of  driftmeter  measurements  is  given  in  Table  2  for  an  8-sec 
period.  The  analysis  algorithm  assumes  that  the  difference  amplifier  is  perfectly 
linear.  The  outputs  of  the  log  electrometer  are  only  output  to  the  TM  once  per 
4  sec  as  LLA  and  LLB  as  8  bit  words  in  SUBCOM.  In  addition,  during  frames 
124  to  127,  LLA  is  output  32  times  per  second  as  a  9  bit  word. 


3.3  Data  Analysis 

The  following  is  an  algorithm  for  converting  the  TM  values  of  the  drift  meas¬ 
urements  (IDM)  into  scientific  units: 

Step  1:  Acquire  a  frame  of  data.  Obtain  the  value  of  COUNT;  this  will 

have  to  be  created  for  odd-numbered  frames.  If  no  more  data  are 
available,  go  to  Step  9. 

Step  2:  If  COUNT  =  N*8,  where  N  =  0,  1,  2 . 15;  Then  to  to  Step  4. 

Step  3:  If  COUNT  =  N*8  +  m,  where  N  =  0,  2,  4 . 14;  Then  to  to  Step  8. 

If  COUNT  =  N*8  +  m,  where  N  =  1,  3,  5 . 15;  Then  go  to  Step  9. 

m  =  1,  2,  3,  ....  7  . 

Step  4:  Calculate  and  retain  the  re-zero  and  offset  values  using  Eqs.  (15) 

and  (17).  The  subscript  denotes  the  IDM  word  within  a  frame. 

*  -  2.52 

(tan  a)H  2  =  2  *  GL  *  k  *  J  *  4**(3  -  Rg) 


(tan  a) 


* 

V,  4 


D4  -  2. 52 

2*GL*k*J*  4**(3  -  R4) 
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asurcment  Word  (I DM) 


pasuroment 


Udii  a;R  g  -  2  *  G.  *  k  *  J  *  4**(3  -  RJ 
L  6 

Dg-2.52 

(tan  a)y>  g  =  2  *  G.  *  k  *J  *  4**(3  -  R.) 

L.  o 

D.  -  2. 52 
2  *  J  *  4**(3  -  Rr> 

Dft  -  2.52 

O  =  - 2 - 

V  2  *  J  *  4**(3  -  Rg)  • 

where 

G.  =1.0  and  k*J  =  -1.  524665. 

L 

If  N  =  1,  3,  5 . 15;  Then  go  to  Step  7. 

N  =  0,  2,  4 . 14.  For  i  =  9  to  16,  calculate  and  retain  the 

drift  measurement  using  Or  from  Step  4  and  Eq.  (21) 

(tan  a)fJ  {  =  (Or  -  (D.  -  2.  52)/(J  *  4**(3  -  R.)))/(GL  *  k) 

Go  to  Step  1. 

For  i  =  9  to  16,  calculate  and  retain  the  drift  measurement  using 
Or  and  Qy  from  Step  4  and  Eq.  (21). 

For  i  =  9,  11,  13,  15 

(tan  a)H  .  =  (Or  -  (D.  -  2.  52)/(J  *  4**(3  -  R.)))/(GL  *  k) 

For  i  =  10,  12,  14,  16 

(tan  a)y  {  =  (Oy  -  (D.  -  2.  52)/(J  *  4**(3  -  Ri)))/(GL  *  k) 

Go  to  Step  1. 

Calculate  and  retain  the  drift  measurement  using  Or  and  Oy  from 
Step  4  and  Fq.  (21) 

(tan  a)R  .  =  <Or  -  (D.  -  2.  52)/(J  *  4**<3  -  R.)))/(GL  *  k) 

Go  to  Step  1. 


Step  9:  Calculate  and  retain  the  drift  measurements  using  Oj.j  and  Oy 

front  Step  4  and  Kq.  (21) 

For  i  =  1,  3,  5 . 15 

(tan  a)jj  .  =  (C>H  -  (D.  -  2.  52)/(J  *  4**(3  -  H.)))/(GL  *  k) 

For  i  =  2,  4,  6,  .  .  . ,  16 

(tan  a)y  .  -  (Ov  -  (D.  -  2.  52)/(J  *  4**(3  -  R.)))/(Gj  *  k) 

Go  to  Step  1. 

Step  10:  The  values  of  the  ion  density  can  be  calculated  as: 

IA  =  10**l(LLA  -  V10)/G,  ) 

IB  =  10**[<LLB  -  V20)/GL] 

Ni  =  (IA  +  +  V'^effJ  * 

Step  11:  There  are  now  four  time  series  of  tangents  of  the  drift  angles  plus 
the  ion  density.  The  tangents  of  the  drift  angles  from  the  standard 
operation  of  the  driftmeter  are  (tan  a)B  and  (tan  a)y.  The  drift 
angle  tangents  obtained  when  the  H+  ions  are  excluded  from  the 
collector  plates  are  (tan  a>B  and  (tan  a)y.  If  satellite  attitude  and 
ephemeris  data  are  available,  the  angles  can  be  converted  to  drift 
velocities  using  Eqs.  (6)  and  (7). 

4.  RETARDING  POTENTIAL  ANALYZER  (RPA)  DATA  PROCESSING 
4.1  Theory  of  Operation 

The  retarding  potential  analyzer  (RPA),  as  shown  in  Figure  4,  is  a  planar 
ion  collector  that  allows  ions  to  freely  enter  the  aperture.  Since  the  spacecraft 
is  supersonic  with  respect  to  all  ion  species,  all  ions  will  pass  through  the 
aperture.  The  ground  plane  around  the  sensor  ensures  that  the  surfaces  near 
the  sensor  aperture  are  electrostatically  uniform,  so  that  straight  lines  represent 
the  ion  trajectories  as  they  enter  the  sensor.  Internally,  th-  retarding  grids  are 
stepped  through  a  series  of  appli'  J  voltages.  When  the  potential  barrier  is  equal 
to  the  energy  gained  by  a  species  from  the  s.n’elliti  not er. tin!  an  1  svelliG-  motion, 


TO  AMPLIFIER 


Figure  4.  Cross-section  of  Retarding  Potential  Analyzer  for  DNA / 
HILAT  Satellite 
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half  the  ions  of  that  species  will  not  penetrate  the  retarding  grids.  For  H  ,  this 
is  0.  29  V  plus  the  satellite  potential.  For  0+,  this  is  4.  59  V  plus  the  satellite 
potential. 

Because  the  satellite  potential  could  go  to  -14  V  with  respect  to  the  plasma 
due  to  exposed  contacts  on  the  solar  panels,  the  ion  sensors  and  ground  plane  have 
been  electrically  separated  from  the  spacecraft  ground  potential.  The  ground 
plane  is  isolated  from  the  sensor  and  electronics  by  a  10^  ohm  resistance.  The 
ground  plane  floats  to  a  potential  with  respect  to  the  plasma  of  0  to  -2  V.  The 
electronics  for  the  sensor  unit  measures  the  potential  between  the  ground  plane 
and  the  spacecraft  (Ab  =  < t>sc  -  -A  floating  potential  of  -A<j>  is  applied,  which 

acts  as  the  ground  potential  for  the  sensors. 

The  current  flow  (I)  to  the  RPA  can  be  expressed  as  a  function  of  the  applied 
potential  (<j> p).  This  analysis  assumes  that  a  single  Maxwellian  temperature 
applies  to  all  ion  species. 


I  - 


AeVo 


T. 


N. 


a.  exp(-xp 


1  +  erf(x.)  + 


(22) 


2 

91  _  -A  e  a 

V2  ffkT 


E 

i  =  1 


(23) 


where 


1  ai 


V  -  ^ (2  eb)/m. 


for  <j>  >  0 


xi  =  V/ai 

X 

erf(x)  =  (2 /N%)  y"  exp(-z^)  dz 


for  $  <  0 


(24) 


(25) 


<j>  =  potential  of  retarding  grid  with  respect  to  plasma  =  b  +  b  (V) 

s  p 

bg  =  potential  of  sensor  aperture  with  respect  to  plasma  (V) 

<j>  =  potential  of  retarding  grid  with  respect  to  aperture  (V) 

2  2 

A  =  sensor  aperture  area  (m  )  =  ^r  ,  where  r  is  aperture  radius  =  0.0127  m 

_1Q 

e  =  electronic  charge  =  1.  602  X  10  C 
a  =  sensor  transparency  =  0.  590 
N.  =  ambient  ion  density  of  i  constituent  (m~  ) 
rm  =  mass  of  i1*1  constituent  (kg) 

=  (2kT’/md1^  =  most  probable  speed  of  i*^  constituent  (m  sec  *) 

T  =  average  ion  temperature  of  all  constituents  (°K) 
k  =  Boltzmann  constant  =  1.38054  X  10  ^  J/°K 
j  =  number  of  ion  constituents 

Vg  =  vehicle's  speed  (m  sec  *)  normal  to  the  sensor  aperture  =  Vg  •  X 

Vj  =  component  of  the  plasma  drift  normal  to  the  sensor  aperture. 

V  =  V  ,  +  V'  . 
d  s 

For  a  spatially  uniform  plasma,  the  current,  I,  is  constant  with  respect  to 
bp  until  x.,  for  the  lightest  ion  present  in  a  measurable  quantity,  approaches 
zero.  Then,  the  current,  I,  decreases  monotonically  and  increases. 

When  is  zero,  3I/9bp  is  at  a  maximum.  This  can  also  be  expressed  by 
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a2i(x.  ^  o)  'a<?2  =  o  .  (2'i) 

i  VP 


becomes  a  constant,  until  x.  for  the  next  heaviest  species  approaches  zero.  If 
there  are  species  present  for  which  |  |  -<1.0  when  x.  =  0,  then  it  will  be  very 

difficult  to  analyze  the  RPA  data.  This  is  the  case  in  the  ionosphere  w.  ->n  both 
H  and  He+  are  present  in  significant  amounts.  If  |  x.  +  ^  |  »  1.  0  when  x.  =  0,  it 
is  relatively  easy  to  find  from  the  data  the  values  of  dp  for  x^  =  0  and  x.+  j  =  0. 
This  is  the  case  in  the  topside  ionosphere  when  and  0+  are  the  only  significant 
ions.  Kxamples  of  HPA  data  are  shown  in  Figure  5  when  Il+  and  0+  are  present 
at  840-km  altitude. 
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Figure  5.  An  Example  of  Ion  RPA  and  Spherical  Electron  Sensor 
Data  as  a  Function  of  Applied  Voltage  at  840-km  Altitude.  The 
upper  plots  are  from  the  spherical  electron  sensor  and  are  typical 
of  Langmuir  probe  data  in  a  large  Debye’length  with  respect  to  the 
size  of  the  instrument.  The  lower  plots  are  from  the  RPA  when 
it  is  in  plasma  composed  of  H+  and  0+ 


In  order  to  find  the  drift  velocity  component  parallel  to  the  sensor  aperture, 
Eq.  (241  is  used.  The  value(s)  of  dp,  where  x.  =  0  is/are  determined  from  the 
data.  If  two  or  more  values  of  dp(x.  -  0)  are  available,  Kq.  (24)  can  be  solved 
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directly  for  the  two  unknowns,  A  and  V  ..  If  only  one  value  of  6  (x.  =  0)  can  be 

s  d  p  \ 

obtained  from  the  RPA  data,  then  A  must  be  obtained  from  another  source,  such 
as  the  electron  sensor. 

If  the  plasma  is  not  uniform  over  the  time  interval  of  a  sequence  of  the 
data  analysis  of  the  RPA  data  will  be  difficult  or  impossible.  Large  variations 
in  the  ambient  density  (up  to  An/n  =  0.  5  at  1G  Hz)  can  be  eliminated  from  the  data 
when  1  Xj|  »  1  for  all  species.  Near  values  of  Xj  =0  in  the  data  variations 
greater  than  An/n  =  0.02  at  16  Hz  can  ruin  the  data  analysis. 

4.2  Operation  of  RPA 

The  aperture  grids  of  the  RPA  are  held  at  the  sensor  ground  potential  at  all 
times.  The  sensor  ground  potential  is  biased  away  from  the  spacecraft  potential, 
as  explained  in  Section  4.  1,  to  keep  it  close  to  the  undisturbed  plasma  potential. 
The  retarding  grids  are  stepped  through  a  sequence  of  applied  voltages  that  are 
shown  schematically  in  Figure  6  and  given  in  Table  3.  The  collector  plate  is 
biased  from  the  sensor  ground  by  the  retar  ding  voltage  plus  -10  V.  The  surpres- 
sor  grid  is  biased  by  the  retarding  voltage  plus  -30  V.  Any  photoelectrons  or 
secondary  electrons  generated  on  the  collector  plate  should  be  driven  back  to  the 
collector  plate  by  the  potential  from  the  surpressor  screen.  Thermal  electrons 
entering  the  aperture  are  prevented  from  reaching  the  collector  by  the  potential 
from  the  suppressor  grid.  Auroral  electrons  generally  have  sufficient  energy  to 
reach  the  collector. 

The  retarding  grid  applied  voltage  is  stepped  through  a  pre-programmed  set 
of  levels  on  a  64-sec  sequence.  At  the  beginning  of  each  sequence,  the  applied 
voltage  is  set  to  the  first  step  in  the  program  (location  00,  0.000  V).  After  each 
readout  of  the  log  electrometer,  the  applied  voltage  is  stepped  to  the  level  given 
in  the  next  program  location.  For  seconds  0  through  3  and  32  through  3  5,  the 
applied  voltage  program  steps  through  locations  00  to  7F.  For  seconds  4  through  5 
and  36  through  37,  the  program  steps  through  locations  80  to  BF  and  repeats 
each  2  sec  (four  frames)  until  second  32,  when  the  first  sequence  is  repented,  or 
until  second  62.  At  second  62,  the  electron  sensor  begins  its  applied  voltage 
sequence.  Design  limitations  have  r.  quired  that  the  RPA  applied  voltages  be  the 
same  as  the  electron  sensor  during  these  two  seconds.  The  voltages  in  locations 
CO  to  FF  are  stepped  through  during  seconds  62  and  63.  At  the  end  of  second  63, 
the  program  re-starts  at  second  0  and  location  00.  If  there  is  a  power  interrup¬ 
tion  at  any  time,  the  program  will  re-start  at  second  0  and  location  00  at  the 
beginning  of  the  next  frame  after  POWKK-ON.  During  testing  it  was  found  that 
the  applied  voltage  program  is  often  one  frame  or'  1/2  sec  ahead  of  the  instrument 
sequencer,  COUNT. 


Table  3.  Applied  Voltages  for  DNA/H1LAT 
e:  Program  locations  are  given  in  hexadecimal) 


Table  3.  Applied  Voltages  for  DNA/IIILAT  (C'ontd) 
Mote:  Program  locations  are  given  in  hexadecimal) 
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c.  t  c.  <r.  c  *r.  c  in  o  ~ 

n-  iM  -H  J.  X  1-  t.3  t  CM  ^  C.  X 


C!©QOOMCOCO»»COI>t-^M,-^r*^ 


ww^nucc^i>^HOcvji>oococoi-T2C'»TOincMt-co’Troc;(>j^'-<,-»7 

OH0)QC0)OOa^OiftC;»OQC5O00©»Oir.iflQC:»OC; 

NOO(Mf-rt*«^a3MO^in^omo^f'^ino^OTt«o^c^ 


--  c:  30  i 


^  Tf  *t  CO  W  M  iM  ' 


I  o  o  o  o  o  < 


KMCMCM'OCO^tiOL'bOwNt-X: 


o^HP5^rtcc-coo<lwUQ^iOHM«^if;ct*coo<;C'J 


CNJCO^i-^t>COO<;KU 


7  JL  OOMncOODOmHCOOLOffiCOCOlfiN^CO^cOOCOOOCOCO^HiOinO 

.2  ~  OOO’t^COOOOOlOOCOOOOOJCOOOOCOC'- 1  o  o  o  oo  c  o  o 

*-*  O^HCvjiOCO^HOnO'^COTfCO'^0(^CNjCn)30C5--*OrHCnCN3I>COTt<r-CNt^CMCO 

C,^  «••■•«•<•••••  •  ••••••  ••  •  • 

Cl  o  ooooO'HHwwwcoco'tin  m*c  ^ocobcM^wHOcocoh-t-ccin 


o-^fMco^in^c-ooo^PQCjQ^ttO^cvico^to^c^ooci^ffiCjQ  —  k* 

COCOCOCOCOOOOOCOCOCOCOWCOCOCOCOOffiOOQOOOQaiCjC-jQOOO 


The  current  to  the  collector  plate  is  measured  by  a  log  electrometer  in  the 
electronics  package.  If  the  current  to  the  plate  becomes  negative  for  any  reason, 
the  output  of  the  electrometer  is  zero.  The  calibration  of  the  electrometer  is 
given  in  Table  4.  There  are  no  modes  of  operation  for  the  log  electrometer.  The 
observed  variations  in  current  collected  are  due  to  the  retarding  potential  and  the 
geophysical  environment.  During  seconds  62  and  63,  the  RPA  and  its  electrometer 
are  functioning,  but  the  data  are  dropped  from  the  TM  output. 

Table  4.  Calibration  of  RPA  Log  Electrometer 


+1 

IQ'12 

o 

1 

»— * 

o 

1 

o 

' 

10'9 

IQ’8 
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1.  26 

2.28 

3.29 

4.  29 

2 

... 
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1.53 

2.  59 

3.  59 

4.  58 

3 

— 

0.  70 

1.74 

2.77 

3.  77 

4.  76 

4 

--- 

0.  82 

1.  86 

2.  89 

3.89 

4.  88 

5 

0.  00 

0.  97 

2.  01 

2.  99 

3.99 

4.98 

6 

0.  02-0. 08 

1. 06 

2.  10 

3. 07 

4.07 

-- 

7 

0.  08-0. 18 

1.  12 

2.  16 

3.  14 

4.  13 

-  - 

8 

0.  12-0.20 

1.  16 

2.20 

3.  19 

4.  19 

-- 

0.  18-0. 26 

1.22 

2.25 

3.24 

4.23 

-- 

0,22-0.30 

1.25 

2.32 

3.29 

4,29 

-- 

!.  .  Current  input  vs  analog  voltage  output  (ambient  temp  =  +14.  8°C  I 

temp  monitor  =  2.  16  V)  1 

. 

+35 

.  1°C  Ambient  temp 

-9.  6°C 

2. 

96  V  Temp  monitor 

1.  14  V 

. 

Input  Current 

Output  Voltage 

Input  Current 

Output  Voltage 

5  X  10'7A 

4.  97  V 

5  X  10'7A 

4.  98  V 

5  X  10'8A 

3.  98  V 

5  X 10'8A 

3.  98  V 

5  X  10"9A 

2.98  V 

5  X  10*9A 

3.00  V 

*' 

5  X  10"10A 

1.91  V 

5  X  10"10A 

2.  01  V 

f 

5  X  10"UA 

0.  92  V 

5  X  10'UA 

1.  02  V 

• 

■ 

-12 

5  X  10 

0.  00  V 

5  X  10*12A 

0.  00  V 

4.3  Data  Analysis 


The  only  accurate  way  to  find  the  drift  velocity  parallel  to  the  sensor  aperturr 
(Vj)  is  to  do  a  least-squares  fit  of  the  data  to  E'q.  (22).  Such  fitting  procedures 
tend  to  require  a  large  amount  of  computing  time.  This  should  not  be  done  in  the 
preliminary  data  processing.  Instead,  only  an  estimate  of  should  be  obtained 
during  the  preliminary  processing  of  the  RPA  data. 

For  an  estimate  of  parameters  of  interest,  it  is  assumed  that  either  «.)+  is 
the  only  species  present  or  H+  and  0+  are  the  only  species  present.  The  follow¬ 
ing  algorithm  is  a  possible  method  of  estimating  results  from  the  RPA  data: 


Step  1:  Obtain  <?g  from  the  Electron  Sensor-  data.  If  unavailable,  assume 
ds  =  -2  V  if  the  vehicle  is  in  the  earth's  shadow,  or  o  ,  =  -0.  5  V  if 
the  vehicle  is  in  sunlight. 


Step  2:  Convert  the  RPA  current  values  (RPA)  from  Tlil  to  log^(l).  Convert 
the  RPA  sweep-voltage  monitor  values  (SWP)  from  TM  to  voltages 
for  frames  0  to  123.  Compare  SWP  values  to  Table  3.  If  the  values 
compare,  match  all  the  currents  with  the  appropriate  applied  poten¬ 
tials  directly  from  the  table.  If  the  values  do  not  compare  exactly, 
adjust  the  applied  potentials  from  the  table  before  using  them.  This 
only  needs  to  be  done  once  per  pass  or  once  per  day. 

Step  3:  Find  maximum  and  minimum  current  values.  If 

(log10(IMAX)  -  log10<lMIM))  is  less  than  5.  then  the  data  are 
inappropriate  for  processing.  Set  all  output  to  some  default  values 
that  indicate  a  null  result.  If  greater  than  1.  5,  use  only  the  data 
points  between  these  two  values. 


Step  4: 


Make  sure  that  log^(I)  is  a  monotonically  decreasing  function  with 
respect  to  increasing  values  of  applied  voltage.  Increases  of  one 
or  two  bits  due  to  random  variations  should  be  smoothed  out. 
Increases  of  five  or  more  bits  from  one  data  point  to  the  next  as  a 
function  of  increasing  applied  potential  should  be  considered  as  too 
great  for  processing.  Set  output  to  default  values  that  indicate  a 
null  result. 


Step  5:  We  will  use  the  fact  that  generally  V  »  a;.  Thus,  the  value  of  the 

expression  in  the  inner  bracket  of  Eq.  (22)  is  ~  2.0  for  A  =0  and 

P 

—  1. 0  for  Xj  =  0, 


First  it  is  necessary  to  find  if  light  ions  are  a  significant  portion 
(>10  percent)  of  the  plasma.  If  significant  light  ions  are  present, 
set  q  to  an  applied  voltage  where  the  light  ions  have  been  completely 


retarded  before  reaching  the  collector,  but  0+  has  not  been  retarded. 
Otherwise,  set  ^  =  0  V, 

Step  6:  Find  bj  such  that. 


logigKbj)  =  iog10I(b0)  -  0.301  . 

Step  7:  From  Eq.  (24) 

mO+(Vs+Vd)2  =  2e(*s+*l)  ‘ 

If  6  is  known,  then  V  ,  can  be  calculated. 
vs  ’  d 

If  bg  is  not  known,  but  the  magnetic  latitude  is  less  than  60°,  then 

we  can  assume  V  ,  ~  0.  A  value  for  6  can  be  calculated  and  used 
d  s 

for  subsequent  data  when  the  magnetic  latitude  is  greater  than  60  . 
If  there  is  a  significant  portion  of  H+  in  the  plasma,  a  second  equa¬ 
tion  like  the  previous  one  can  be  found  by  approximately  the  same 
technique.  Then,  the  two  equations  can  be  used  to  simultaneously 

solve  for  V  ,  and  b  . 

d  s 

Step  8:  Find  <l>^ (>bj)  such  that. 


ioglO1^  *  logjo1^!)  -  0.319 
From  Eq.  (24) 


2kT/mQ,  =  (  2 


V  -  yj(2e<t>9)/m  ^ 

^  o+J 


5.  ELECTRON  SENSOR  DATA  PROCESSING 
5.1  Theory  of  Operation 

The  two  ion  sensors  (Driftmeter  and  RPA)  are  variations  of  a  Faraday  Cup 
design.  This  means  that  all  components  of  the  plasma  pass  through  the  aperture 
with  as  little  change  in  their  ambient  trajectories  as  possible.  Once  inside  the 
sensor,  the  components  are  separated  by  electrostatic  potentials.  The  electron 
sensor  is  a  Langmuir  probe,  which  is  characterized  by  an  electrostatic  potential 
being  exposed  to  the  ambient  environment  and  is  capable  of  drawing  current  from 
the  environment.  Under  this  loose  definition,  the  exposed  interconnections  on  the 
solar  panels  are  Langmuir  probes.  A  scientifically-useful  Langmuir  probe 
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measures  the  current  to  a  collector  as  a  function  of  the  electrostatic  potential  on 
the  collector  in  a  controlled  manner. 

There  are  two  basic  types  of  Langmuir  probes,  a  cylindrial  probe  and  a 
spherical  probe.  The  Electron  Sensor  is  a  modified,  single  spherical  probe. 

A  spherical  probe  is  spherical  in  shape  and  it  is  assumed  that  the  sheath  around 
the  probe  is  spherical.  The  sheath  is  the  depth  into  the  plasma  that  an  electro¬ 
static  potential  penetrates.  Since  an  electrostatic  potential  falls  with  a  charac¬ 
teristic  depth  of  the  Debye  length  (X^)  from  the  source,  the  sheath  is  an  arbitrary 
number  of  Debye  lengths  (typically  3  *  to  10  *  X^).  For  830 -km  altitude,  the 
Debye  length  can  be  expected  to  range  from  0.2  to  11  cm.  To  ensure  that  the 
sheath  is  spherical,  the  probe  should  be  removed  from  all  other  parts  of  the  space¬ 
craft  by  the  sheath  depth.  At  high  densities,  this  will  be  true.  At  low  densities, 
there  may  be  some  problem  with  the  sensor  sheath  being  affected  by  the  spacecraft 
surfaces. 

The  two  types  of  spherical  probes  are  the  double  floating  probe  and  the  single 
floating  probe.  In  the  double -floating  probe  system,  two  probes  are  electrostat¬ 
ically  biased  with  respect  to  each  other.  The  current  passes  from  the  plasma  to 
one  probe,  through  the  measuring  system  to  the  second  probe  and  back  to  the 
plasma.  A  single  probe  system  replaces  the  second  probe  with  the  reference 
potential  or  ''ground1’.  The  spacecraft  single  probe  uses  the  sensor  ground  poten¬ 
tial,  which  is  related  to  the  spacecraft  ground  by  a  floating  bias  potential.  Since 
the  spacecraft  ground  is  not  capable  of  supplying  infinite  current,  the  Electron 
Sensor  can  affect  the  spacecraft  potential  if  the  sensor  attempts  to  draw  too  much 
negative  current. 

The  Electron  Sensor  has  two  spherical  elements.  The  outer  spherical  ele¬ 
ment  is  a  spherical  grid  to  which  a  potential  is  applied.  The  inner  element  is 
a  solid  sphere  to  which  a  potential  of  <j>  +  20V  is  applied.  The  current  to  the 

inner  sphere  is  measured.  This  configuration  assures  that  thermal  ions  are  not 
a  component  of  the  collected  current  and  it  minimizes  the  effect  of  photoelectron 
emissions  from  the  collector. 

The  negative  current  (-1^)  to  the  spherical  collector  is  a  function  of  the  poten¬ 
tial  of  the  sensor  grid  with  respect  to  the  plasma  <#.  If  <t>  is  negative,  electrons 
that  approach  the  sensor  are  retarded.  If  <j>  is  positive,  electrons  within  the  sheath 
are  accelerated  toward  the  sensor.  If  </>  is  zero,  the  electron  current  to  the  sen¬ 
sor  is  the  "random  current".  Assuming  a  single  population  of  electrons  with  a 
Maxwellian  distribution  function,  the  current  of  electrons  passing  through  the 
spherical  grid  is  given  below: 
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for  tfr  <  0, 

I  =  -A  e  a  a  exp  (-x2)/(2\/tr)  (27) 

for  <f>  >  0  and  »  r, 

I  =  -AeoN  a(l  +  x2)/2^  ,  <28> 

e  e 

where 

x2  =  e</>/kTe 

<j>  =  potential  of  sensor  grid  with  respect  to  plasma  =  4>s  +  (V) 

<j>g  =  potential  of  sensor  ground  with  respect  to  plasma  (V) 

$  =  potential  of  sensor  grid  with  respect  to  sensor  ground  (V) 

^  2  2 
A  =  surface  area  of  sensor  =  3.24  if  r  =4.  104E-4  (m  ) 

r  =  radius  of  the  spherical  grid  =  6.  35E-3  (m) 

a  =  sensor  grid  transparency  =  0.  648 

a  =  (2kTe/m  )*/2  =  most  probable  electron  speed  (m  sec"1) 

Tp  =  electron  temperature  (°K). 

Ideally,  the  results  from  the  Electron  Sensor  swept  applied  voltage  sequence 

would  follow  the  solid  curve  in  Figure  7.  To  the  left  of  <j>  =  -<f>  ,  Eq,  (27)  would 

P  ^ 

describe  the  data.  Much  of  the  data  will  approximai  ely  fit  these  equations,  but 
not  all  of  the  data.  Notice  that  in  the  retarding  region,  Eq.  127)  results  in  a 
straight  line  in  semi-log  coordinates 

a  log10(-Ip)/a«p  =  5040°K/Te  .  (29) 

There  are  several  reasons  that  the  measured  current  may  not  exactly  match 
the  theoretical  current.  First,  there  may  be  a  hot  component  to  the  electron 
population  as  well  as  a  cold  component.  This  is  most  likely  to  occur  in  the  sub- 
auroral  region  where  the  local  ionosphere  is  dark  and  the  conjugate  ionosphere  is 
sunlit.  The  total  current  is  the  sum  of  two  components,  each  of  which  can  be 
fitted  by  Eq.  (27).  In  such  a  case,  T^  =  1500°K,  T2  *  25,  000°K,  and 

Nj  2  103  *  N2. 

Another  source  of  departure  of  the  measured  current  from  the  theoretical  is 
locally-produced  electrons  that  reach  the  collector.  Figure  7  shows  the  measure¬ 
ment  of  secondary  electrons  caused  by  accelerated  ions  (mostly  He  )  striking  the 
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grid  wire.  In  sunlight  there  will  be  a  small  current  du 

on  the  grid  wire  and  sensor  Doom  reaehing  the  coileeto; 

The  most  troublesome  source  of  1, pasture  of  the  rn 

theoretical  current  is  when  X  and-'or  c  change  while  p 

O  '  s 

can  occur  due  to  spatial  gradients  in  the  ionospheric  pla 
getic  particles  impacting  the  spacecraft. 


to  photoelecl ton 


ensured  current 
is  being  varied 
sma  and/or  due  ■ 


APPLIED  POTENTIAL  (VOLTS) 


Figure  7.  Example  of  log^Q<-Ip)  vs  $  for  HILAT 
presents  theoretical  re 


Electron 

Sensor.  Solid  line  represents  theoretical  response.  Dots 
represent  data  points  during  a  sweeping  applied  voltage 
sequence 
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5.2  Operation  of  Electron  Sensor 


The  output  from  the  Electron  Sensor  consists  of  the  output  from  the  log  elec¬ 
trometer  and  from  four  passband  filters.  As  shown  in  Figure  6,  the  Electron 
Sensor's  grid  is  held  at  a  potential  of  +1.  5  V  with  respect  to  the  sensor  ground  for 
62  sec  of  each  64-sec  period.  During  the  first  62  sec,  the  output  ELEC  represents 

_  n 

a  dc  sampling  of  the  current  collected  from  approximately  6  X  10  amps  to 
6  X  10  ^  amp.  During  the  63rd  and  64th  seconds,  the  electron  sensor  has  a  series 
of  potentials  applied  from  +4.  725  V  to  -4.  725  V  based  on  locations  from  CO  to  FF 
in  the  stored  program  given  in  Table  3.  The  calibration  of  the  Electron  Sensor- 
log  electrometer  is  given  in  Table  5. 

The  four  passband  filter  outputs  represent  the  ac  output  of  the  Electron  Sensor 
log  electrometer  after  passing  through  a  series  of  passbands.  The  four  filter  out¬ 
put  are  centered  at  70,  220,  70,  and  2200  Hz.  Figure  8  shows  the  bandwidth  of 
each  filter.  The  output  of  each  filter  is  proportional  to  the  rms  of  Alog^I-I^). 
Assuming  that  is  constant,  the  filter  output  is  also  related  to  rms  (AN/N). 

Table  6  shows  the  relationship  between  the  output  signal  and  Alog^f-I^  ).  During 
the  63rd  and  64th  seconds,  the  outputs  from  the  filters  are  not  meaningful  due  to 
the  applied  voltages. 

5.3  Data  Analysis 

5.3.  1  DENSITY  MODE  ANALYSIS 

The  majority  of  Electron  Sensor  data  is  taken  when  the  applied  potential  is 
static.  We  will  assume  that  the  potential  with  respect  to  the  plasma  is  either 
accelerating  or  zero.  Thus,  variations  in  current  to  the  Electron  Sensor  are 
assumed  to  be  due  only  to  geophysical  density  variations.  The  measured  electron 
density  is  then 

log10Ne  *  log10(’V  '  logio(A  *  e  *  °  *  aA2  *  •  (30) 

For  density  fluctuations  in  the  range  16  to  0  Hz  can  be  found  by  taking  the  rms  of 
log10^e  at  t*le  aPPr0Pr’at;e  interval: 

Alo*10(‘V  =  rmS  tlog10('Ie)i  -  log10(-Ie)i+j)  •  J  =  1,  2,  3,  .  .  . 

o-  log10(ANp/Np)  .  (31) 
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INPUT:  de.=  IC37  A. 


AMBIENT  TEMP*»  14.8* C 
TEMP.  MONITOR  *2.16  V 


a  a=  4876  <07Cl  ( F I ) 

28.l5mVrms*l07fl(F2) 


Input:  d.C  *  10  7A 

a  e.*  16-25  mVrms  *  I0u  (F3) 

9-386  mVrms.  l07fl(F4) 


FREQUENCY  (Htl 

Figure  8.  DN'A/HILAT  Langmuir  Probe  Filter  Bank  Response  Functions 
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should  be  scaled  so  that  N^  =  N.  when  both  value  are  available.  If  the  sealing 
factor  is  greater  than  10,  the  Electron  Sensor  data  should  be  considered  as  invalid 
due  to  vehicle  potential  problems  (d  <  -2.  5  V). 

5.3.2  SWEPT  APPLIED  POTENTIAL  MODE  ANALYSIS 

The  analysis  of  the  Electron  Sensor  data  during  seconds  62  and  63  as  a  func¬ 
tion  of  applied  voltage  yields  the  sensor  ground  potential  (dg),  the  electron  tem¬ 
perature  (T  )  and  absolute  electron  density  (N  ).  There  is  no  correction  of  the 
e  '  e 

applied  voltage  to  the  voltage  in  the  grid  spacings. 

The  following  algorithm  has  proved  successful  in  the  past  and  is  recommended 
for  the  present  data  set.  It  is  not  an  exact  match  of  the  theoretical  formulae 
(Eqs.  (27)  and  (28)).  but  it  is  quick  on  a  digital  computer  and  it  considers  the  dis¬ 
crete  digitization  of  the  data.  The  algorithm  is  as  follows: 

Step  1:  The  applied  voltage  monitor  output  (SWP)  is  checked  again  the 

expected  values  in  Table  3.  If  the  values  check,  the  applied  voltage 
values  from  Table  3  with  a  -4.  725  V  bias  are  matched  with  the 
measurements  of  electron  current  (ELEC).  If  the  values  do  not 
match,  an  adjustment  is  made  in  the  stored  values  to  match  the 
actual  values  and  then  they  are  paired  with  the  measured  electron 
current. 

Step  2:  The  measured  values  of  log^l-I^)  vs  dp  are  searched  for  "glitches" 

(obvious  random  noise  in  the  data  stream)  that  should  be  ignored. 

Step  3:  The  maximum  and  minimum  currents  are  found.  All  data  not 

between  <f> p ( I )  and  tfp(I^jj^)  are  discarded.  All  currents  in  the 
saturation  level  of  the  instrument  are  discarded. 

Step  4:  If  the  first  four  points  after  the  maximum  current  are  fitted  to  a 

straight  line,  this  forms  the  dashed  line  in  the  acceleration  region 
of  Figure  7. 

Step  5:  The  remaining  data  are  searched  to  find  the  best  straight  line  for 

the  retarding  regions.  This  search  is  done  by  making  a  straight 
line  fit  through  each  group  of  four  points  to  find  the  maximum 
absolute  value  of  the  slope  of  log,„(-I  )  vs  d  .  If  two  or  more  fits 
give  equal  maximum  slopes,  take  the  value  for  the  most  positive 
value  of  bp-  Let  the  slope  of  this  line  be  "S". 

St  -p  6:  Find  the  intersection  of  the  straight  lines  from  steps  4  and  5.  The 

potential  at  the  intersection  (d  °)  is  the  zero  potential  with  respect 
to  the  plasma,  which  in  turn  gives  the  potential  of  the  sensor- 
ground: 
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The  potential  of  the  vehicle  with  respect  to  the  plasma  is 


</>sc  =  os  -  SEN POT  .  (33) 

Step  7:  The  electron  temperature  is 

Te  =  5040°K/S  .  (3  4) 

Step  8:  The  absolute  electron  density  is 

N  =  I  (<j>  °)  *  2  *  v'tr / (A  *  e  *a  (2  *  k  *  T  /m  )1'2)  .  (3  5) 

6  0  p  1  0'  e 

6.  CONCLUSIONS 

The  preliminary  data-processing  output  for  the  Thermal  Plasma  Experiment 
should  consist  of  three  elements: 

(1)  Raw  data  without  unpacking  or  conversion  from  TM  units  to  scientific 
units.  These  data  must  be  time  tagged  correlation  with  other  experiments  and 
with  attitude/ephemeris  data. 

(2)  Raw  data  unpacked  and  converted  to  scientific  units. 

(3)  Preliminary  results  calculated  from  the  data. 

The  information  for  unpacking  the  raw  data  is  given  in  Section  ?.  The  con¬ 
versions  from  TM  units  to  scientific  units  are  summarized  in  Table  7. 

The  calculations  of  preliminary  results  are  summarized  in  Table  8. 


Table  7.  Summary  of  Conversions  of  Raw -Data  From  TM  Units 


Word 

Reference  for  Conversion 

COUNT 

None 

RPA 

Table  4 

SWP 

Table  3 

IDM 

Section  3.3  and  Table  2 

LLA,  LLB 

Eq.  (9)  and  Section  3.3 

ELEC 

Table  5 

TEMP 

Eq.  <1> 

SENPOT 

Eq.  (2) 

FI,  F2,  F3,  F4 

Figure  8,  Table  6,  and  Eq.  (32) 

Table  8.  Summary  of  Preliminary  Calculated  Results 


Value 

Samples  per  sec 

Frames  Available 

Reference 

N.(DM) 

1 

0-123 

Section  3.  3 

l 

32 

124-127 

N 

32 

0-123 

Eq.  (30) 

e 

1/2 

124-127 

Eq.  (35) 

3/2 

8-63,  72-123 

Eqs.  (22) 

H  O 

1/3 

0-5,  64-69 

to  (26) 

1 

6-7,  70-71 

AN  / N 
e'  e 

1 

0-123 

Eq.  (32) 

70,  220,  700,  2200  Hz 

ANe/Ne 

— 

0-123 

Eq.  (31) 

(16/j)Hz,  j  =  l,2,  3... 

V,,  V 

16 

all 

Section  3.  3 

~h'  ~v 

*s 

1/64 

124-127 

Eq.  (33) 

Vd(DM) 

16 

all 

Eq.  (6) 

Vd(RPA) 

... 

... 

Eq.  (24) 

List  of  Symbols 


Scale  variables 
Vector  variables 
Unit  vector 

Components  of  unit  vector 
Multiplication  of  scalars 
Scalar  multiplication  of  vectors 
Exponentials 


A  or  a 
B  or  x 

X 

Bx  or  x 
A*B  or  2a 
A  •  B 

“v*  *>- 

A**2  or  A2 


